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The propagation of hard partons through spatiaUy extended matter leads 
to medium-modifications of their fragmentation pattern. Here, we review 
the current status of calculations of the corresponding medium-induced 
gluon radiation, and how this radiation affects hadronic observables at 
collider energies. 
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1. Introduction 

In the QCD improved parton model, inclusive cross sections for higli-pj^ 
partons (observed as leading hadrons or jets) can be calculated through 
coUinear factorization. This formalism convolutes the perturbatively cal- 
culable parton-parton cross sections with non-perturbative but process- 
independent ("initial state") parton distribution functions and ("final 
state") fragmentation functions. What additional physics input is needed if 
one extends the application of QCD to the description of inclusive high-pj^ 
parton cross sections in nucleus-nucleus collisions at RHIC or LHC, where 
hadronic projectiles have significant spatial extension oc A^/"^ ? 

The current discussion of this question focuses in particular on the fol- 
lowing effects: (i) The density of the incoming parton distributions increases 
for given x and by a factor ^ .X^i^. im plies that with increasing 

A^^^ the non-linear modificationsEmanEnESIof the QCD evolution equa- 
tions become relevant at lower center of mass energies (or larger momentum 
fractions x). (ii) The QCD evolution of the projectile wave function may be 
altered due to the presence of a spatially extended "target" through which 
it propagateJ^. This goes under the name "initial state partonic energy 
loss" . It is discussed e.g. in Drell-Yan production in p — A where a quark 
in the proton may un dergo multiple scattering in the nucleus on its route 
to annihilation (jji) xhe hard partonic collision is well- localized on a 

scale 1/Q which is much smaller than the diameter of a nucleon 1/Aqcd- 
Thus, one does not expect medium-modifications of the hard partonic cross 
section itself. However, there is a class of high-pj^ observables in e-A and p-A 
for which the medium-modification can be described by nuclear enhanced 
twist-four matrix elements whose importa nce is de termined by the pole 
structure of the hard partonic cross section B^ | 44 | 5U | ^j^jg jg ^j^g only case 
of a nuclear dependence for which higher twist factorization theorems are 
proven, iv) At sufficiently large center of mass energy or for sufficiently ex- 
tended projectiles {A^^^ > 1), the probability of two independent hard par- 
ton interactions in the same hadronic collision becomes significant^!. This 
is e.g. a significant background to some of the hadronic Higgs decay chan- 
nels searched for in p-p collisions at the LHC (v) Finally, the presence of 
spatially extended matter can affect the fragmentatioiiandhadronization 
of hard partons produced in nucleus-nucleus collision J 1 '^ l '^4 | 3| xhis effect 
is often referred to as "final state partonic energy loss" or "jet quenching" 
and is expected to affect essentially all high-p_L hadronic observables in 
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heavy ion collisions at collider energies. The main motivation for the study 
of these observables is that the degree of nuclear modification may allow 
for a detailed characterization of the hot and dense matter produced in the 
collision. 

All the above mentioned effects should emerge from a theory of QCD 
processes in nuclear matter. The development of such a theory is still at 
the very beginning. The high rate at which new data are becoming avail- 
able from RHIC, as well as the rapid development of new theoretical ideas 
and phenomenological interpretations make it likely that any comprehen- 
sive review will be outdated very soon. In this situation, we have decided 
to focus entirely on a simple partonic multiple scattering formalism which 
underlies at present almost all "jet quenching" calculations. This limited 
scope excludes many important issues. In particular, we do not discuss to 
what extent ex peri mental data from RHIC provide indications of partonic 
energy loss (seelsSl). Moreover, many phenomenologically proposals for dis- 
covering p artonic en ergy loss, such as the measurement of back-to-back 
correlationfl3SESl the specific behavior of massive &-quarks ^ 42 i l7l 
or particle ratios at high pj_ are not included in the discussion. Finally, we 
do not discuss recent alternative formulations of partonic energy lo ss, e.g. 
in the framework of nuclear enhanced higher twist matrix elements l^v l vv l 

The presentation is organized as follows: in section|21 we discuss the for- 
malism of eikonal scattering in which 5-matrix amplitudes are determined 
in terms of eikonal Wilson lines in the target field. Section |31 discusses the 
generalization of this formalism to non-eikonal trajectories which take the 
transverse Brownian motion of the partonic projectile into account. The 
medium-induced gluon energy distribution radiated off a hard parton is 
calculated in this generalized formalism. Section ^discusses the properties 
of this medium-induced spectrum and sectionl^lturns to some applications. 



2. Gluon Bremsstrahlung in the Eikonal Approximation 

At high energy, the propagation time of a parton through a target is short, 
partons propagate independently of each other and their transverse posi- 
tions do not change during the propagation. The only effect of the propa- 
gation is that the wave function of each parton in the proie ctile a cquires 
an eikonal phase due to the interaction with the target fielc^^^. To be 
specific, we consider a hadronic projectile wave function whose relevant de- 
grees of freedom are the transverse positions and color states of the partons, 
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^-zn^ V'({a»,xJ) |{a,,x,;}) . (1) 

{Qi:Xi} 

The color index can belong to the fundamental, anti fundamental or 
adjoint representation of the color SU (N) group, corresponding to quark, 
anti quark or gluon in the wave function. We choose the Lorentz frame in 
which the projectile moves in the negative z direction. Thus at high energy 
the light cone coordinate x+ of the projectile does not change during the 
propagation thought the target. The projectile emerges form the interaction 
region with the wave function 

^out^^^'™^ V'({a.,x,})[]VF(x,)a.ft|{A,xJ). (2) 

{ai,Xi} i 

Here S is the 5'-matrix, and the W^s are Wilson lines along the (straight 
line) trajectories of the propagating particles 

WiK,) = V exp{i J dz-T"A+(x„ z-)} (3) 

with - the gauge field in the target and T" - the generator of SU{N) in 
a representation corresponding to a given parton'^. 

The interaction in the target field changes the relative phases between 
components of the wave function and thus "decoheres" the initial state. As 
a result the final state is different from the initial one, and contains emitted 
gluons. 

The component of the outgoing wave function which belongs to the 
subspace orthogonal to the incoming state reads 

\6^) = [1- |*™)(^',;„|]|*o„t). (4) 

From this, one can calculate for example the inelastic cross section which 
up to the factor of total flux is given by the probability 

Pincl = = 1 - |SP , (5) 

^In this section we use the light cone gauge A~ = 0, in which the parton distributions of 
the projectile are simply expressed in terms of the parton number operators. In this gauge 
the color fields of the target have a large component, thus the eikonal S - matrix 
eq.J2l is given in terms of A+. The physical observables are of course gauge invariant 
and can be calculated in any gauge. In Appendix A we present the same calculation in 
the "target light cone gauge", A+ = 0, where the only non vanishing component of the 
target field is A^. 
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where 

S=(*m|*o«t>. (6) 

The S'-matrix element in eq.® has to be averaged over the distribution 
of the color fields in the target. Also, the total and diffractive cross sections 
can be calculated by projecting other components out of '^out 

In what follows, we apply the above framework to calculate the gluon 
radiation spectrum radiated off a hard quark which propagates at high 
energy through a nuclear target. At this point we assume the quark to be 
coming from outside the target rather than being produced inside the target 
in a hard scattering event. Thus it carries with it the fully developed wave 
function which contains the cloud of quasi real gluons. 

In the first order in perturbation theory the incoming wave function 
contains the Fock state \a) of the bare quark, supplemented by the coherent 
state of quasi real gluons which build up the Weizsacker- Williams field /(x), 

^ l«) +/ dxdC/(x)r„%|/3;fe(x,C)) (7) 

,b 



u a , g ^ „ 



Tap 

Here Lorentz and spin indices are suppressed. In the projectile light cone 
gauge A~ = 0, the gluon field of the projectile is the Weizsacker- Williams 
field 

AXx)cx0(x-)/,(x), /,(x)«5^, (8) 

x^ 

where x~ = is the light cone coordinate of the quark in the wave function. 
The integration over the rapidity of the gluon in the wave function |(7J) goes 
over the gluon rapidities smaller than that of the quark. In the leading 
logarithmic order the wave function does not depend on rapidity and we 
suppress the rapidity label in the following. 

The interaction of the projectile (0) with the target leads to a color 
rotation /3j of each projectile component i, resulting in an eikonal 

phase W{yii)aifji- The outgoing wave function reads 

*out = W[.,{Q) |7) + j d^fi^) T^pW^^m W,ti^) |7 ; c(x)) , (9) 



Gluon Radiation and Parton Energy Loss 



7 



where T^^(O) and M^^(x) are the Wilson hnes in the fundamental and ad- 
joint representations respectively, corresponding to the propagating quark 
at the transverse position Xg = and gluon at Xg = x. Gluon radiation is 
an inelastic process, for which we have to calculate the projection 
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(10) 



W^p(0) 



a 



Wpy(O) 



a 



Here, the index 7 in the projection operator runs over the quark color index. 
It appears in ea- lfTU)! since we have to project out the components of the 
wave function which differ from the initial state only by a global rotation 
in the color space EHt. 

The number spectrum of produced gluons is obtained from H10|l by cal- 
culating the expectation value of the number operator in the state S'^a, 
averaged over the incoming color index a. After some color algebra, one 
obtains 

^prod(k) - 4 V(<5vl/„|4(k) a,(k)| 5^^) 



27r 



dx dy e 



ik.(x-y) X ■ y 
y 
1 



■jj-2 y2 



1 



^^_-{(Yr[W^\y)W^{Q)])), 



1 



-{{Tr[W^\y)W^{^)]))^ 



(11) 



Here, we have used /(x) /(y) — ^ '^yi for the Weizsacker-Wilhams field 
of the quark projectile in configuration space and the symbol ()t denotes 
the averaging over the gluon fields of the target. 

To obtain an explicit expression, we require the target average of the 
product of adjoint Wilson lines. Various averaging procedures are discussed 



''An event in which the final state difli'ers from the initial one only by a global rotation 
of the color index would be an elastic scattering event. 
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in the literature l O-'^ l -'-^ l 'J < | og | jo | specific, we consider a target which 

consists of static scattering centers with scattering potentials (q) at po- 
sitions (x„,z„). In the high energy approximation, each scattering center 
transfers only transverse momentum to the projectile, 

At{^,z-)^Y. J (^e^(™>'ia^(q)5(z--z-). (12) 

We define the target average as an average over the transverse positions of 
the static scattering centers. Introducing the longitudinal density of scat- 
tering centers, n{z^) = S{z^ — z~), we find 

{{J dz-dz-A+i^,z-)A+{y,S-)))t= j din{i)^a{^-y), (13) 

where we have defined the dipole cross section 

a(x-y) = 2| (0|a+(q)P {l-e^^-^^y^) . (14) 

We now consider the gluon number spectrum (|ll(l in two limiting cases. 



2.1. N—1 opacity approximation (single hard scattering 
limit ) 

This limit is relevant for a small target with weak target fields v4+ — 0{g). 
In this case the projectile can scatter at most once. Thus the adjoint Wilson 
lines can be expanded to leading order in the single scattering potential 
a+(q). Using we find 

-l—{{TT[W^H^)W^{y)]))t^l-^ Jd^niOai^-y). (15) 
The corresponding gluon number spectrum reads 

f7(x) +a{y) -CT(x-y) 

Differentiated with respect to the momentum transfer q, this result is the 
well-known Gunion-Bertsch radiation cross section^ for the gluon radia- 



A^prod(k) 



8tt 



d^ n(^) dxdy e 



»k.(x-y) x-y 
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tion cross section in quark-quark scattering, 



da 



GB 



Ca as 



+ 0(1/E) 



(q)l' 



(17) 



dlnxdqdk 7r2 '"^^iJ^-" k2 (k - q)2 • 

Integrated over the transverse momentum exchanged with the "target" 
quark, this radiation spectrum can be written in terms of the dipole cross 
section H14II . 



d In a; dk 



= AC A as 



(2 7r)2 k2 (k-q)2 



-1 



dr 



— I (^:;;^^QCD(r) exp{iq-r} 



(18) 



For a scattering on a static source described by a Yukawa type potential. 



one has ' 



*QCD 



(q)l^ 



Although perturbatively gluons are mass- 



less, the screening mass M is frequently introduced to mimic the medium 



effects which cut off the infrared divergence^. The factor as in this expres- 
sion accounts for the fact that the gluon field in eas. p3ll4|) is emitted by 
the single particle source. It is the presence of this coupling constant that 
justifies the expansion of the Wilson lines to leading order in a. 

It is interesting to compare the expression eq. I|17|) with the high energy 
limit of the Bethe-Heitler photon radiation cross section in QED, 

)(q)l 1,2 /I, ^„S2 (19) 



d\nx dq dk tt-^ 



'QED I 



k2 (k - xqy 



k,co 



'^The absolute normalization of a+ in eQS. 1171 and 1161 is different, which accounts for 
the different constant factors in these equations. 
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Here, k is the transverse photon momentum and x = lu^/E the fraction 
of the energy carried away by the photon. Again, if integrated over the 
transverse momentum transfer q from the target, this radiation spectrum 
can be written in terms of a dipole cross section 



The radiation cross section IjlTI) and (|19|l for QCD and QED differ only 
by their x-dependence. The QCD radiation is flat in rapidity, whereas the 
QED radiation cross section peaks at projectile rapidity x — 1. The physical 
reason is that in QCD, the radiated gluon is charged. The distribution of 
gluons in the incoming quark wave function is flat in rapidity. Thus when 
a gluon scatters directly off a scattering center into the final state, the 
rapidity plateau is produced. It is the non-abelian diagram contributing 
to (|17|l which dominates the QCD radiation spectrum. This is also the 
reason why the explicit Casimir factor in (|17() is the adjoint and not the 
fundamental one. In QED on the other hand only the electron can scatter 
directly off a source. This scattering is most effective for electrons with 
small fraction of longitudinal momentum. Such scattering events result in 
final state photons which carry most of the initial longitudinal momentum, 
and thus the photon spectrum is peaked around x = 1. 

In configuration space, the different x-dependence of the QCD and QED 
radiation spectrum is reflected in the different arguments of the dipole cross 
sections: the QCD dipole cross section is tested at transverse distances r but 
the QED one is tested at a;r. For an intuitive understanding of this behavior, 
consider the incoming projectile electron or quark in the light cone frame 
as a superposition of the bare particle and its higher Fock states, 



If all Fock components interact with the scattering potentials with the same 
amplitude, then the coherence between these amplitudes is not disturbed, 
and no bremsstrahlung is generated. The radiation amplitude depends on 
the decoherence between the different Fock state elements, i.e., on the dif- 
ference between the elastic scattering amplitudes of different fluctuations. 
This difference is characterized by the impact parameter difference between 
different Fock state components. It can be estimated from the transverse 
separation of the \q) and |g7) fluctuations. On a characteristic time scale 




(20) 



Iprojectilc) = \q) + \qg) + ... . 



(21) 
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'Ke ~ (^ 2x{i-x)E ^ ^^'^ transverse energy of the two-particle Fock 

state, the quark, gluon or photon move in the impact parameter plane by 
distances 

fc_L 1 (l-x) 
x^,,.:^,oc— — oc^— , (22) 

Since the photon does not carry charge, the transverse distance x±-y plays 
no role in the decoherence of QED Fock state components, while the dis- 
tance x± g plays the dominant role in the corresponding QCD process. The 
QED and QCD dipole cross sections thus measure the transverse separa- 
tions 

1^ (l-x)r ■►^-xr--^| 

QED: Y e e 

QCD: g q q 

To leading order 0{l/x), the cross sections in ()18(l and H2()|l thus differ by 
a factor x in the argument of cr(r). Based on the above picture, Nikolaev 
and Zakharov proposecpSl to write the dipole cross section for the qqg- 
configuration beyond leading 0{\/x) as 

t7QCD(r,x) = ^ {a(r)+a((l-x)r)}~ia(xr), (24) 

where the prefactors account for the color representations of quark and 
gluon. 

Due to this distinct x-dependence the QED expression cannot be ob- 
tained in the eikonal approximation. The eikonal limit where the energy 
loss of the electron is neglected is equivalent to the limit a; = 0, where 
according to ea. H19() no photons are emitted. In our wave function formal- 
ism this arises since the components of the wave function with no photons 
and with one photon rotate with the same phase during the propagation 
through the target, if the energies of the electron in both components are 
assumed to be the same. Thus in QED the eikonal approximation leads to 
the vanishing inelastic cross section. In QCD this is of course not the case, 
since the gluon component of the wave function scatters most efficiently. 



12 



A. Kovner and U.A. Wiedemann 



2.2. Multiple soft scattering limit (dipole approximation) 

We now consider the opposite limit, when the projectile undergoes many 
scatterings. Assume that these scatterings are independent in the sense 
that the vector potentials of the different scattering centers in eq. H12I) are 
uncorrelated in color space. Then the dipole cross section in ea. H15|) expo- 
nentiates. We also assume that the cross section is dominated by the low 
momentum transfer in each individual scattering. In this limit the dipole 
cross section (t(x — y) can be approximated by the lowest order term in 
the Taylor expansion a{x — y) cx (x — y)^. This is the so-called logarith- 
mic approximation, in which the target correlator of two Wilson lines is 
characterized in terms of the saturation momentum Qs, 



Ar2 



exp 



exp 



J d^niOai^-y) 

Ql 



Ca 

8 Cf' 



(25) 



Here, the saturation momentum is 

-,2 f dq 



Qi 



i2ny 



I4(q: 



|2 „2 



(26) 



which in the model H12|) of the target fields characterizes the average trans- 
verse momentum squared transferred from the target to the projectile. Us- 
ing the av erag e (^5)1 to calculate the produced gluon number spectrum l(TT|l. 
we obtain 

271 j x^y2 



1 



(2^)3 Ca Ql ^ 



dqe 



CaQ 



k2 (k - q)2 



(27) 



Expression H27|) contains again the Gunion-Bertsch radiation cross section 
(|17|l . In the dipole approximation (|25|l . however, the entire medium acts 
coherently as a single scattering center with a Gaussian momentum distri- 
bution 



cx e 



(28) 
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There is an interesting attempt to generalize (127 |l to the case of a dense 
projectile realized in nucleus-nucleus collisions'"^. 

The eikonal approximation described in this section is valid as long as 
the average momentum transfer from the medium is much smaller than 
the longitudinal momentum of the produced gluon. Only in this limit all 
scattering centers contribute coherently to the production amplitude. The 
purpose of the following is to generalize H27|) to the kinematic region where 
this is not necessarily the case, and incoherent scatterings also play an 
important role. 



3. Gluon Radiation beyond the Eikonal Approximation 

In the eikonal approximation, the Wilson lines Q determine the ^-matrix 
for the scattering of partonic wave functions. In this section, we derive 
the leading finite energy corrections 0{l/p~) to this high energy limit. 
We explain the form of the target averages including 0{l/p~) corrections 
and we give two applications: the photo absorption cross section and the 
medium-induced gluon radiation spectrum. 



3.1. Wilson Lines for Non-Eikonal Trajectories 

We consider the A^-fold scattering amplitude of a high energy quark. In the 
light cone gauge A" = 0, it reads 

F(a;o,p) = 



p ,a p ,a p ,a p ,a 

»> .-^ .-^ 



P .a 



N=0 



N 



0© 

dpf dp~ dpj 



(27r)4 



dxt dx- dx, 



— m + le 



X (-Z7- A+(x,,a;-)) e"^(P^+l-f> )+«r-(pii-p+)e»x,.(p,+i-p.)^ 

xe"°Piw(p). (29) 

Here, the target gauge fields = A+T^ are path-ordered along the tra- 
jectory of the quark; we suppress color and spin indices on the quark prop- 
agators and the spinor v{p). The additional phase factor 6"" ^^ in H29|l is 
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useful in order to connect the multiple scattering diagram to the vertex 
of a more complicated amplitude by integration over xo- In the absence of 
scattering centers, (|^ reduces to a free quark wave function e^'^°'Pv{p). 

We approximate eq. to leading order O in the norm and 

to subleading order O (l/p^) in the phase. To this order we can neglect the 
dependence of the target fields on xf , as indicated in H29|l . This is because 
the incoming quark propagates in the negative z direction and thus only 
feels the fields at x^ — 0. The integration over xf and p~ is then trivial: all 
propagators have the same component = P - To simplify H29|) further, 



^sUs{p)uM ^2_^ (30) 



^ — m P ITT- 2p 

and neglect the last term to leading order in the energy l/p^ . The spinor 
structure can be simplified: 

Us{p) (-«7 • A) Us'ip) ~ A+ Sss> ■ (31) 
With these simplifications, eq. 1)29(1 takes the form 

oo N 



- (p? + "i^)/p 

xe'^(!'-P"e""P^"*^''-Pu(p). (32) 



The spatial ordering of the scattering centers along the trajectory is fixed 
explicitly by performing the p^-integrations. Neglecting the mass term, we 
find a product of 0-functions times phase factors 



dpf i 
2tt p+ - (pf + m2)/p- + ie 



e(xr - xT_^) . (33) 



Inserting this expression into H32I) and performing the transverse p^- 
integration, we obtain the basic building block of the formalism, the free 
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light-cone Green's function 



(27r)2 



-ipi ■ (xi_ 1 — Xf 



I?r(^) exp 



exp 



(34) 



Here, r = dv/dz^ and the boundary conditions on the free path integral 
are v{x~) — x^. The limit x'^ — x^_^ — > of H34|l is equivalent to the high- 
energy limit 



lim Go(xj_i, a:-_;^;xi,a;- \p ) ^ 5 (xj_i - 



(35) 



The free Green's function H34|l evolves the free plane wave e'^^ ^ -«xi P 
the position {yii-i,x~_^). The sum of the iV-fold scattering amplitudes 
reads 

oo / N „ 

F{xo,p) = e"o P" H n / c'2;,^dx,Go(xi_i,a;:ii;x„a;7|p-) 

JV=0 \ 1=1 



xe(xr-xr^) (z.4+(x„x-)) e^iv^^" 



-IXjv-p 



u(p) 



z— ^ — zpi -x 

(ixG(xo,a;Q ;x,x"|p")e w(p) , (36) 



where the full Green's function is 
G{xo,Xq;x, x^\p^) 

T^HO exp ^ 
W{[r];xQ ,x^) = 7^ exp 



d^r\0 



W{[v];x^,x-), (37) 



z / d^A+{r{0,0 



(38) 



These expressions give the leading 0{l/p^) corrections to the phase of 
the eikonal Wilson line The path integral in 1)3 7|l has a very simple 
interpretation. It describes the quantum mechanical particle which, while 
propagating through the medium can move in the transverse plane. The free 
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transverse motion is governed by the kinetic term - a two-dimensional 
light-cone hamiltonian with "mass" given by the total energy p~ /2 of the 
projectile along the beam. During the motion the wave function of the 
particle also acquires a phase given by the Wilson line along the trajectory. 
The infinite energy limit corresponds to the classical limit, where a single 
straight classical trajectory dominates the sum over paths. In this limit the 
Green's function (|37|l reduces to the eikonal Wilson line in 

lim G(xo,Xo ;r,x"|p") = Ty(r;xo ,0:") 5^(xo - r) . (39) 

3.1.1. Non-abelian Furry approximation (target field Aq) 

The multiple scattering diagram contains the leading order 0{l/p^) 
corrections to the phase of the eikonal Wilson line, but the norm of (|36() 
is given to O ((l/p^)°) only. For a scalar quark that would be the only 
important correction. For a spinor quark however this is not enough if one 
wants for example to discuss emission of on shell gluons. It turns out that 
the leading order of this propagator vanishes when connected to an emission 
vertex of a real gluon. To get a leading order result for this process one has 
to keep order 0{l/p~) corrections to the first propagator 1/(7 • pi — m). 
We now discuss the non Abelian Furry approximation which addresses this 
issue. 

We also change our notations slightly in accordance with literature. 
The target field will not be described b y a larg e A^-component but rather 
specified by the Gyulassy-Wang model l ^^ l ^^v l as a collection of scattering 
potentials of Yukawa-type with Debye screening mass A/, 

A^(x)=<5o^^(p^(x)^^ (40) 

i=l 

^^(x) =^(x-x,)^"°-- (41) 

Here, the i-th scattering center is located at x^ and exchanges a specific 
colour charge a — Oi. 

"^Intuitively one thinks of large A+ as being a more appropriate description of a target 
which itself is moving in the positive X3 direction, while large being more appropriate 
for a static target. Technically, however there is no difference between the two sets of 
vector potentials as long as the projectile is very energetic and its x+ coordinate does 
not change during the propagation inside the target. In this sense switching to large Aq 
can be considered as a different gauge fixing. 
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For this target field (|40|l , the calculation of the multiple scattering dia- 
gram H29(l proceeds in close analogy to section mi The components 7^^+ 
in are replaced by 70^0- The role of the p^-integrals is played by 
the ^/-integrals which lead to an ordering along the longitudinal direction. 
To simplify the spinor structure, one uses (^)7o — £■(70 ~ 73)70 which can 
be iterated, (£(70 - 73)70)" v{p) = 2""! E'^ (70 - 73)70 v{p) ~ 2" v{p). 

In comparison to section 13.11 the only new feature of the present dis- 
cussion is that we approximate the first quark propagator to 0{1/E). In a 
coordinate system with p || n and for o? = 70 7, we have 



^70 = p{lo - 73)70 + a-^ • (pi - p) - ^ ^ Je 



p(7o - 73)70 + ia 



dy 



a ■ (p — pnj 



(42) 



Here, the operator ^ acts on wave functions of the form exp{— ip; • y}. 
This structure leads to the differential operator Di 



V a-(p, -npj) 



2E 



2E 



(43) 



which carries the transverse and longitudinal spinor structure to order 
0{1/E). To lowest non-trivial order 0{1/E) in the norm and phase, the 
multiple scattering diagram takes now the for 



^'«(2/°,y,p) = e 



DF{y,p)u{p), 



(44) 



r 1' I 



P ,a 

M 3 



P ,a 



N=0 



O © © 



P ,a 

MSI N 



p. a 



Here, the differential operator D acts on the transverse wavefunctions F 
which for far forward longitudinal distances xl satisfies plane wave bound- 
ary conditions 



Foo{-^,XL,Pt) = exp <^ -ipi • x + i —'-xl > , (45) 
F{y,yL,Pz)^ / dxG(y,?/L;x,a;L|p)Foo(x,a;L,p) . (46) 
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The full Green's function in this expression is given in close analogy to 
(|87|l by the sum over paths weighted with the Wilson lines. It satisfies the 
recursion relation involving the free non-interacting Green's function Go, 

G{r, z; r', z'\p) = Go(r, z; r', z'\p) 



-i I di I dpGo{T, z] p,£,\p) Ao{p,0 Go{p,i]v' , z'\p) 



XL XL 



+'P j dil j d^2 J dpidp2Go{r,z;pi,^i\p)iAo{pi,^i) 
X Go (pi , Ci ; P2 , 6 b) « ^0 (P2 , 6 ) G(p2 , 6 ; r' , z» , (47) 



Go(r,z;r ,0 \p) = 



■ exp ■ 



ip (r — r') 



(48) 



2m{z' - z) 1 2(z' - z) 

The path ordering V in (|47|l ensures that the potentials Ao are ordered 
along the longitudinal axis. 

This formula is in fact a direct non Abelian generalization of the cor- 
responding expression in Quantum Electrodynamics. Consider the solution 
of the Dirac equation for an electron in a spatially extended scattering 
potential C/(x) 

d 



I ^ — U{yi) — m 7o + i a • V 



*(a;,p) ^ 0. 



(49) 



In the high-energy approximation in which we neglect corrections of order 
0{U/Ei) a.nd 0(1/Ef), the solution to this Dirac equation is the Furry 
wave function I I \ 

^f{x, p^) = e^^-*-»IP-l^ A F{k, p,) v{p,) . (50) 

This solution has the same form as the non-abelian multiple scattering 
diagram in H44|l . The Green's function used to evolve the wave function F 
is now given in terms of the scattering potential U{x), 

G(r,z;r',z'b) = 



VriO exp <^ / 



|r2(0-*C/(r(0,0 



(51) 



Again, it describes the quantum mechanical evolution of the projectile par- 
ticle in the plane transverse to the beam. 
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3.2. Target averages for non-eikonal Wilson lines 

In section 121 wc discussed the target average of two eikonal Wilson lines. In 
section IXTl we showed that the 0(l/i?)-improved generalization of eikonal 
Wilson lines are given in terms of the path integral H37I) . We can compute 
the target average of two such Green's functions expanding both Green's 
functions in powers of Aq according to (|47() and then resumming the series. 
We assume that the different scattering centers in (|40|) are uncorrelated 
and that the projectile can interact with a given scattering center only once 
(elastically or inelastically). For generality we take the two Greens' function 
to act on initial states with different energies (1 — a)p and ap, respectively. 
This is a generic situation in the application to deeply inelastic scattering, 
where the initial quark and anti quark share the energy in the initial state. 
This will also be useful for us in the following, since we will need to calculate 
phase differences between the wave functions of an initial state parton and 
a final state parton with a different energy. Performing the ta rget average 
as in H13(l . the target average for two adjoint Green's functions'^ leads to 
the appearance of the dipolc cross section ° a, 



((G(ri,z;r'i,z'(l - a)p)|G'(r2, z; r'j, z'|ap)»t 
Vri / I?r2 exp 



ri(z)=ri 



r2(2:)=r2 
1 



de|(«r2 + (l-«)r?) 



+i-n(0cT(ri(0-r2(e)) 



(52) 



With the help of the transformation 



r,(0 = (l-a)ri(0 + ar2(0, 



(53) 
(54) 



the target average can be written as 



((G(ri, z; r'l, z'(l - a)p)|G(r2, z; r'„ z'\ap)))t 

= /Co(ra(z'), z'; ra{z),z\p) JC{rb{z), z'; rb{z),z\pa{l ~ a)) . (55) 



"^In the rest of this article we use cr to denote the cross section for scattering of dipole 
consisting of two adjoint charges. It differs from a defined in by the factor 
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Here, the new path integrals are of the form 



/Co(r',z';r, ^Im) 



2tt i (z' — z) 



exp 



i/i (r' — r) 



2(z' - z) 



Z 

/C(r', 2'; r, z|/i) — J X>r exp i J 



|r2 + zin(0a(r) 



(56) 
(57) 



In the high-energy hmit, this path-integral coincides with the average (|26|) 
of two eikonal (fundamental) Wilson lines, 



lim /C(r', z'; r, z|/i) = exp 

fl — *oo 



-n(Oa(r) 



5(2)(r'-r) 



(58) 



The path- integral H57|) determines the 0(l/i?)-corrections to this eikonal 
average. Quite generally, it is seen to compare the relative distance (|54|) 
between the trajectories of the two generalized Wilson lines as a function 
ofe 



3.3. The medium-induced gluon radiation spectrum 

We consider now the medium-induced radiation amplitude for a quark of 
color a splitting into a quark and a gluon of colors b and c respectively. The 
diagram for arbitrarily many M-, L-, and TV-fold scatterings of the in- and 
outgoing partons is 



® © 



© 



© © 



© © 



© 



(59) 



These multiple scattering contributions are summed up in the Green's func- 
tions which include the 0{1/E) corrections to the eikonal Wilson lines. 
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The corresponding radiation amplitude reads 



Ma^bc ^-ij dyi e*«^' e"^!^'! j d^y j d^^i d^^2 d^iCg 

2 2 2 

^gixl■pl-^X2■p2-^Xg■kg-^^x^+^^x^-^^x^ 

where the interaction vertex reads 



(60) 



(61) 



Here, the differential operators Di stem from the Furry approximations of 
the quark lines, see eq. H43|l . The gluon radiation amplitude (|60|l is a typical 
example for which the order 0(i?°) term in the norm of the Green's function 
vanishes. As a consequence, the derivatives ^ which are subleading in the 
Green's function itself, appear to leading order in the transition amplitude 
and in the radiation cross section. 



1 



d(lnx)dk (27r)2 AEl{l-xy 



dp {\Ma- 



(27r)2 x^Ef{l-xy J^_ Jy^ 
X J dydy dxi dxi dug d5ig dx2 e^^^'^"^^^"^^^ . 

x^G^-(xi;ybOT-,^ (|^G-^(y; x,|^)) G^J)^(y; x^b^) 

G^:5(x,;yb2) (|G-5(x,;y|.)) J-^^^G^(j;5c,\p^^ 
Graphically, this can be visualized as 



(62) 



(63) 



As discussed earlier, the softer parton, which in this case is the gluon, 
rescatters much more efficiently. Keeping only the leading 0{x) terms the 
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final result for this radiation cross section readJ^ 

X /due"»'' "e-^-fs<^'^«"(«'^(") ^ . i_/C(y = 0,y,;u,^zM. (64) 
J dy du 

The e-regularization in this expression has been introduced to suppress the 
contributions of gluons emitted at i = ±oo. With this adiabatic switching 
on (off) of the interaction at asymptotically early (late) times, the state 
of a single quark is an eigenstate of the Hamiltonian at t —^ ±oo, and 
thus the medium induced emission of an extra gluon is indeed a transition 
to a different eigenstate. One could in principle perform the calculation 
without introduction of this explicit regulator. However in this case one 
must take care to calculate the transition amplitude between eigenstates of 
the interacting Hamiltonian and thus take into account 0{as) corrections 
to the wave functions of the initial and final states. The initial "dressed 
quark" state in such a calculation will contain the admixture of a \qg) 
component as discussed in Section 2. To calculate the radiation spectrum 
one then needs to subtract from the final state wave function the component 
which overlaps with the "dressed quark" wave function. Introduction of the 
regulator makes the calculation technically much simpler. However, care is 
needed since the limit e — > in (|64() does not commute with the longitudinal 
integrations. In practice, one can remove the regularization e.g. by splitting 
up the longitudinal y^- and j/i-integrations of the radiation cross sections 
into six parts, 

2+ z+ L 2+ L L L 2+ 2+2 + 

(65) 



2_ yL z_ j/_L z_ z_ L yL L L yL 

and to write the radiation cross section as a sum of the corresponding six 
contributions 

^^(^ = ^C.(/i+/. + /3+/4 + /5+/a). (66) 

This allows to remove in all terms Ij the e-regularization of (|64|l . In what 
follows, we are interested in the gluon energy distribution of a quark pro- 
duced in the medium at position z- = 0. In this case, the contributions 
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-^1 5 ^2 J -^3 to (|65|l vanish. The remaining three terms denote the direct pro- 
duction term I4 for gluon emission inside the target in both amphtude and 
complex conjugate amplitude, the destructive interference term which is 
the characteristic interference term between gluon emission in and beyond 
the target, and the hard vacuum radiation term /g. 



14= 



15= 2 Re 



.0 



l6 = 



(67) 

The contribution Iq can be associated to the medium-independent contribu- 
tio n to the hard radiation off a nascent quark jet propagating in the vacuum 
This is a consequence of the leading 0{x) approximation in which only 
the rescattering of the radiated gluon contributes, and the rescattering of 
the hard quark in H()7|l is neglected, 



^6- j^- 



(68) 



This term may be viewed as the medium-unaffected DGLAP part of the 
gluon radiation. In the following, we shall discuss the medium-induced mod- 
ification of the gluon radiation spectrum, determined by and I5. For 
notational convenience, we shall not denote explicitly that the medium- 
independent term /g is subtracted. 



4. Properties of the medium-induced gluon radiation 
spectrum 

Two approximation schemes for the radiation cross sect ion H64|l have 
been explore d rece ntly: the dipole approximation Ej^SEHl and the opac- 
ity expansiorl '^^ l '^^ l In the dipole approximation, the medium-dependence 
results from the multiple soft scattering of the projectile in the spatially 
extended matter. The opacity expansion can be related to an expansion in 
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the effective number of scatterings and allows for hard momentum transfers 
from the medium. Here, we discuss for both approximations the inclusive 
energy distribution of gluon radiation off an in-medium produced parton. 
Integrating over t he tra nsverse momentum of the radiated gluon up to 
|k| < x^^, we haveEHEOI 



dl 
] — 

du! 



d_ 
dy 



2Re / dyi I dy, 



d_ 
du 



2?r exp 



d\ 



d'ke 



-i/~d€«(0<T(u) 



niO<j{r) 



(69) 



The radiation off hard quarks or gluons differs by the Casimir factor Cr = 
Cp or Ca, respectively. The properties of the medium enter eq. by 
the product of the density n{^) of scattering centers times the dipole cross 
section cr(r) which determines the strength of a single elastic scattering, see 
(|14|l . Note that the density can be time-dependent as is for example the 
case in expanding medium. 



4.1. Multiple soft scattering 

Consider a medium-dependence introduced by arbitrary many soft scatter- 
ing centers, rather than a few hard ones. In such a medium the projectile 
performs a Brownian motion in transverse plane; the transverse position 
of the projectile in configuration space changes rather smoothly and small 
relative distances in the dipole cross section are important. The integrand 
of (|64|l has its main support at small transverse distances r — | rj whi ch 
allows to write the dipole cross section to logarithmic accuracy as^^^ 

n(0a(r)^ig(0r^ (70) 

Here, q{^) is the transport coefficienlPl which characterizes the transverse 
momentum squared jjp' transferred to the projectile per mean free path A. 
For a static medium, it is time-independent, 

(n, 

In the multiple soft scattering approximation, the transport coefficient q and 
the in-medium pathlength L are the only informations about the properties 
of the medium which enter the gluon radiation spectrum H69(l . 



Gluon Radiation and Parton Energy Loss 



25 



4.1.1. The harmonic oscillator approximation (static medium) 

In the so-called dipole approximation {TUJ, the path integral (|57ll in the 
gluon r adiatio n cross section (|64|) is equivalent to that of a harmonic 
oscillatoj™ 

j[ 

ICosc{y, yi-u,yi\fi) ^ —exp[iAB{y^ + u'^) - 2i Ay u] , (72) 

2 sm(S2 Ay) 
with complex oscillator frequency 

n^^xl^. (74) 

Explicit solution for parton produced inside the medium: 

Inserting /Cose into (^^, the medium-induced gluon energy distribution 

becomes 

^i^^^Cn{h + h) , (75) 

/j = / difik / fc_L dk± Ij , (76) 
Jo Jo 

where the terms Ij are defined by eq. Ht)7|) . Their explicit form is 
/4 — 47rRe / dzi I dz; 



2 1 - cosh[(l + i)Az] 



4- 1 e^^*^=^ sinh[(l+i)Az/2]/2 I ^yy-j 



/s = 47rRe / dz 



4 

1 + z 



2 sinh[(l + i)5] 

(1 -f i) e^(^+*)*^=^ tanh[(l-(-4)i/2] 



(78) 



4 sinh[(l + i)z/2] cosh[(l -h i)z/2] / ' 
F = 2l~- 2icoth[(l + i)Az/2] [{I + i){Ml + Z2) 

-1- coth[(l -I- i)A5/2]] , (79) 

N = 22 sinh f ^(5i - Z2)) + (1 - z) cosh f ^^^(ii - ^2)) ■ (80) 
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Here, the longitudinal distances are written in dimensionless units, rescaled 
by the modulus of the oscillator frequency |ri|, 

z = V2\n\z, L = V2\n\L. (81) 

The upper integration limit of the k± integral takes in dimensionless vari- 
ables the from 

M, = ^ . (82) 



'2ujV2\n\ 

Totally coherent limit for parton in initial state: 

Assume that a parton is part of an initial state. In this case, prior to scat- 
tering on a target, the parton has already a fully developped wavefunction 
with additional virtual gluons. In the scattering process, these virtual glu- 
ons can be freed. As a consequence, the gluon energy distribution radiated 
off this parton has to include radiation vertices at arbitrary early times. 
One needs to include all six terms Ij of Ht)t)|) in calculating eq. H75|l . Ex- 
plicit expressions for this case are given in Ref. EOl. Here, we consider only 
the high-energy limit of the resulting expression. Since fl cx this high- 
energy limit of the radiation spectrum corresponds to the leading order in 
0(L = V^ML). It reads 

dtodk 7r2 qL J (27r)2^ ' (q - k)2 k2 ' ^ ' 

This is the totally coherent limit of the gluon radiation spectrum for an 
incoming parton. It is the Gunion-Bertsch gluon radiation spectrum 128|) for 
a target momentum transfer q which is acquired by the transverse Brownian 
motion of the incoming projectile. 



4.1.2. Qualitative estimates oftu^ vs. quantitative calculations 

Important properties of the gluon energy distribution can be obtained from 
simple estimates. Their range of applicability, however, depends on phase 
space constraints which are not as easy to e stiniat e. Here, we compare 
qualitative arguments to numerical calculationJ^^^ES 

Qualitative argumentsl^ We consider a gluon in the hard parton wave 
function. This gluon is emitted due to multiple scattering if it picks up 
sufficient transverse momentum to decohere from the partonic projectile. 



Gluon Radiation and Parton Energy Loss 



27 



For this, the average phase tp accumulated by the gluon should be of order 
one, 




(84) 



Thus, for a hard parton traversing a finite pathlenth L in the medium, the 
scale of the radiated energy distribution is set by the "characteristic gluon 
frequency" 



(85) 



For an estimate of the shape of the energy distribution, we consider the 
number A'coh of scattering centers which add coherently in the gluon phase 
(I84|l . fc^ ~ A'coh M^- Based on expressions for the coherence time of the 
emitted gluon, Uoh - -p- - a/^ and iVcoh 



for the gluon energy spectrum per unit pathlength 



dl 



1 



tcoh 




one estimates 



duj dz A'coh du) dz icoh " V t^-" ^ ^ 

This 1 / v^-energy dependence of the medium-induced non-abelian gluon 
energy spectrum is expected for sufficiently small lo < lOc- 

Quantitative analysis: The gluon energy distribution l|t)9|l dep ends not 
only on cjc, but also on the dimensionless "density parameter'i^ 



(87) 



A finite value of R implements the kinematical constraint fcj^ < 
the transverse momentum phase space of the emitted gluon, see (|69|l . The 
kinematical boundary is % = 1. This constraint is neglected in the argument 
leading to the l/^/oJ-energy dependence of 

The limit R oo removes the kinematical constraint from Ht)9|) and 
coincides with the original result of Baier, Dokshitzer, Mueller, Peigne and 
Schifli", 



lim cj — 

ij— >oo dui 



2asCR 



In 



The limiting case for small u is as expected from the estimates (|86|l , ^ 
dl 2asCR { ^/W. for 



lim w — 



2a) 



for 



< Wc , 



(89) 




Fig. 1. (a) Dependence of the medium-induced radiative energy loss (AE) on the in- 
medium pathlength L for E = 100 GeV. The transport coefficient is defin ed here as 
^ = no C (b) Double logarithmic presentation of (a). Figure taken fromlHOJ 

The average parton energy loss is the zeroth moment of this energy distri- 
bution 

{AE)r^^ = Um fdcoLO^^ iu,ocL\ (90) 

R-^oo J duj 2, 

This is the well-known i^-dependencc of the average energy lossEEEl. The 
numerical ca lcul ation in Figure ^ shows deviations from the limiting case 
(|90|l . [In an erroneous prefactor Nc was included in the definition 

of cj^; hence, the values for AE shown in Fig^ and 13 are a factor 3 too 
large.] In general, the strong increase of the average energy loss Ai? with 
in-medium pathlength indicates the strong sensitivity of partonic energy 
loss to the geometry of the nuclear collision. 

The average energy loss /S.E is the first moment of the energy distri- 
bution . Fig. 121 shows this distribution for finite values of the density 
parameter R. At sufficiently large gluon energy, the distribution approaches 
for any value of R the BDMPS limit IjSSI) . Below a critical gluon energy 
ill, however, the finite size gluon spectrum is depleted in comparison to 
the BDMPS limit. This suppression results from a pha se sp ace constraint: 
gluons are radiated on average at a characteristic angle' 

k\ V^jq f ^ \ 1 



^<i- (91) 



This depletes the gluon distribution for gluon energy lu below 

LOc \R, 



I \ 2/3 



(92) 
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\ BDMPS (R = oo) 




Fig. 2. The energy distribution of radiated gluons for different values of the kine- 

matical constraint R = lUc L. Figure taken fromE3 

This i?-dependent suppression is seen clearly in Fig. |21 



4.1.3. Angular Dependence of the gluon energy distribution 

Quantum Chromodynamics is a finite resolution theory; the energy of a 
single quark or gluon is not a measurable quantity. One rather measures the 
fraction of the hadronic decay products of a single quark or gluon which are 
radiated within a certain jet opening angle. Thus, the angular dependence 
of the medium-induced modification of gluon bremsstrahlung is important. 

The angular dependence of can be studied by varying the kinematic 
cut-off X of the transverse momentum integrationl ^ l ^^ l^. For fixed charac- 
teristic gluon energy lUc = ^Q^^ and fixed dimensionless parameter ^qL^, 
the radiation spectrum oj-^ inside a given angle = arcsin(x) depends 
on the parameters lOc and R = x^igL^. From Fig. |21one thus concludes 
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immediately that the harder gluons are emitted more coUinear. 

From the gluon energy distribution inside a given cone opening angle, 
one obtains the medium-induced additional average energy emitted outside 
a cone of given opening angle 6, 

A TP A F A 

— (6 = arcsin(x)) = — (x = 1) - -^(x) • (93) 
As seen from Fig. O ^^(B) does not decrease monotonously with in- 




Fig. 3. The fraction of the total radiative energy loss AE/E emitted outside a jet cone 
of fi xed angle ©. Here, the transport coefficient is defined as q = 2noC. Figure taken 

frominni 

creasing Q but has a maximum at finite jet opening angle. The reason is 
that the radiative energy loss outside a cone angle 9 receives additional 
contributions from the Brownian fcj^ -broadening 

k2 ' (k + q)2 ^^^^ 

of the hard vacuum radiation term /g. Such contributions do not affect the 
total k-integrated yield -^(O = 0), since they result only in a shifting of 
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the transverse momentum phase space distribution of the emitted gluon. 
However, this shift in transverse phase space shows up as soon as a finite 
cone size is chosen. We conclude from Fig. El that the total k-integrated 
radiative energy loss^(e = 0) is not the upper bound for the radiative 
energy loss outside a finite jet cone angle -^(O). Depending on the trans- 
port coefficient and the in-medium path length, the latter can be larger by 
more than a factor 2. 



4.1.4. Harmonic oscillator approximation (expanding medium) 

In nucleus-nucleus collisions at collider energies, the produced hard partons 
propagate through a rapidly expanding medium. The density of scattering 
centers and thus the transport coefficient q(^) is expected to reach a max- 
imal value qd around the plasma formation time and then decreases 
rapidly due to the strong longitudinal and - to a lesser extent - transverse 
expansion, 

9(0 = 9^(1)". (95) 

Here, a = characterizes the static medium discussed above. The value 
a = 1 corresponds to a one-dimensional, boost-invariant longitudinal ex- 
pansion and approximates the findings of hydrodynamical simulations. The 
formation time of the medium may be set by the inverse of the saturation 
scale psat and is « 0.2 fm/c at RHIC and « 0.1 fm/c at LHC. Since 
the time difference between the formation of the hard parton and the for- 
mation of the medium bulk is irrelevant for the evaluation of the radiation 
spectrum (|69|) . one can replace in (|69|l the production time of the parton 
by 0. 

For a dynamically evolving medium of the type l|95|l . the path-integral 
(|57|l in H()9I) is the path integral of a 2-dimensional harmonic oscillator 
with time dependent (imaginary) frequency fi^(0 = and mass uj. 
The explicit solution can be written in terms of the variables z(^) = 

2zj.f](^o)^o (£f'^'^ and v ^ 

K.{r{yi),yi;r{yi),yi\uj) ^ (96) 
^ {cir{yiy + C2v{jjif - 2r{yi) ■ v{yi)) 



exp 



2'KiD{yi,yi 



2iD{yi,yi) 
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Fig. 4. LHS: The medium-induced gluon energy distribution radiation for a dynamically 
expanding collision regions I95i with expansion parameter a = 0, 0.5, 1.0 and 1.5. The 
value of the transport coefficient qo is taken at initial time ^o- RHS: The same gluon 
radiation spectrum with parameters rescaled according to HIOOI . Figure taken from ^21 



Here, the coefRcients are given by modified Bessel functions I^{z) and K^{z) 
[here: z = z{yi), z= z{yL)] 

ci = ci(yL,2/L) = z{z/zY [/,_i(z)i^,(z) +i^,_i(z)/,(z)] , (97) 
C2 = C2(yL,2/L) = z{z/zy [K,{z)I,^i{z) + I,{z)K,^i{z)] . (98) 

and the fluctuation determinant reads 

^(^'^") = Tr-^rhr^^'^or [iA^^)KAz^) - k,{z)i,{z^)] . (99) 

One checks that for a = 0, = |, the solution \T1\ is regained. 

Using the exphcit form of the path integral (|96|l . one can calculate the 
medium-induced gluon energy distribution 169(1 for a dynamically expand- 
ing medium ^. The result is shown in Fig^] The radiation spectrum 
satisfies a simple scaling law which relates the radiation spectrum of a dy- 
namically expanding collisio n re gion to an equivalent static scenario. The 
linearly weighed line integrally 

l=jjifj^'^d£,{^-^o)m (100) 
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defines the transport coefficient of the equivalent static scenario. The linear 
weight in HlOOfl implies that scattering centers which are further separated 
from the production point of the hard parton are more effective in leading 
to partonic energy loss. In contrast to earlier believe that parton energy loss 
is most sensitive to the hottest and densest initial stage of the collision, this 
implies for a dynamical expansion following Bjorken scaling [a = 1 in eq. 
(|95|l ] that all timescales contribute equally to the average transport coeffi- 
cient. This makes partonic energy loss a valuable tool for the measurement 
of the quark-gluon plasma lifetime. 

4.1.5. Properties of the transport coefficient 

The only property of the medium which enters the medium-induced energy 
distribution H69() is the transport coefficient 

9 = A^VA = I ^q2|a+(q)p (101) 

For phenomenological estimates of the value q for cold nuclear matter, 
one often invokes the relation of q to the gluon structure function q = 
^^^^^^f:^p[xG{x,qL)]. This leads to estimates S gcoid < (200MeV)Vfm. 

There are alternative approaches EE which extract ^coW from the medium 
modification of hard processes in "cold nuclear matter" , such as medium- 
modifications t o the D rell-Yan production D in p — A or the dijet imbal- 
ance in e — A | 4^ | 44 | r^j^^ extracted values are consistent with (jcoU < 
(200 MeV)^/fm, but the systematic uncertainties of these determinations 
are difficult to specify. 

For the transport coefhcient in hot quark-gluon matter, even less is 
known. Baier, Dokshitzer, Mueller and Schiff 1^ estimated for the rescat- 
tering properties of hot matter at a temperature T = 250 MeV a value 
q — (500 MeV)^/fm. In any case, the transport coefficient is expected to 
increase linearly with the energy density, as seen in Fig. [S] In the case of a 
phase transition at Tc « 175 MeV, this indicates a strong jump of q at Tc- 
The medium modification of hard processes observed at RHIC and LHC 
should allow to determine the value of q and thus provide an independent 
estimate of the energy density attained in the collision. For this, the proper 
treatment of the dynamical expansion of the collision region e.g. via 1)100(1 
is obviously important. 
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Fig. 5. Es tim ate of the transport coefficient as a function of the energy density e. Figure 
taken from^Sl 

The transport coefficient fias been related also to the saturation mo- 
mentum scaled 

Ql^qL. (102) 

The saturation scale entering here should be defined for a color octet dipole 
- thus, it is a factor Ca/Cf larger than the saturation scale determined 
from the qq scattering probability in DIS. The numerical value of Qs is 
very uncertain but < (3 GeV)^ may be considered as an upper bound 
at LHC. For in-medium path lengths L up to twice the nuclear size, this 
corresponds to values up to i? = ^qL^ < 40000. 

4.2. Opacity Expansion of the radiation cross section 

In what follows, we discuss the opacity expansion of the radiation cross 
section H6H|I . This is based on the expansion of the path- integral H57|l in 
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powers of the dipole cross section, 

Z 

-J d^n{C) y"rfri/Co(r,y/;ri,0 ^cr(ri)/Co(ri,^;r,yi) 

Z 

Z Z 

+ j diin{ii) J d^2n{^2) jdvi drj /Co(r, y;; ri, ^i) 

xicr(ri)/C(ri,$i;r2,6) ^cr(r2)/Co(r2,6;f,y;)- (103) 

The iV-th order in opacity is obtained by expanding the int egrand of the 
gluon energy distribution (EH) to A^-th order in {n{£) <7{r))^ | 59 | 26 | 2 5|^ 
start with a discussion of the first order term for which the medium acts ef- 
fectively as a single hard momentum transfer positioned within path length 
L. 



4.2.1. Expansion to order N — and N — 1 

For the gluon energy distribution H69(l of a free incoming quark, the zeroth 
order in opacity vanishes. A free incoming particle does not radiate without 
interaction. A parton produced in some hard process, however, does radi- 
ate without further interaction in order to decrease its virtuality. This is 
the DGLAP parton shower whose medium modification is accessed by the 
medium-induced spectrum lu^. To zeroth order in Ln^, this spectrum is 
given by 



d^I{N = 0) 
duj dk 



®- 



as ^ 1 



(104) 
(105) 



This is the characteristic gluon radiation spectrum associated to the produc- 
tion of a hard parton. Medium modifications to 1)104(1 appear to first order 
in opacity. Introducing the transverse energies Q = j^, Qi — ^^^^r~^ ^'^^ 
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finds 
(Pl{N 



duj dk 



7r2 2w2 
-k - qi 



■ no ■ 



(2^|ao(qi)| 
LQi - sin(L(3i) 



(By 



O 




(E>- 



o 



(106) 



(107) 



The formation time for a medium-induced gluon is tf « l/Qi- For a finite 
distance ~ L between the production point of the parton and the point of 
its rescattering, there are thus gluons whose formation times are too long: 
L/tf = LQi < 1. The interference term ^Q^^^MLQi) prevents 
such gluons from being radiated. It provides a logarithmic cut-off for an 
otherwise infrared divergent k- integrated expression, see eq. (|117|) below. 

This radiation spectrum interpolates between the totally coherent and 
totally incoherent limiting cases. These can be discussed in terms of the iso- 
lated contributions for the hard (vacuum) and the medium-induced Gunion- 
Bertsch radiation terms: 

H(k,)= ■ 



H(k^+ q. 



o 



o 



o 



6 



(108) 



In the totally coherent limit, the extension of the target is negligible, and 
medium-corrections are absent, 



(fa^^^^HN = 1) 
hm --. r-:: = . 



(109) 



L^o d{liix)dk 

In the incoherent L — > oo limit, the radiation cross section 1)69(1 expanded 
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up to first order in opacity, takes the form 




d^I{m) 



doj dk 



-^CF{l-wi)H{k) 



H — ^ Cf no L 



TT- 




The three terms on tlie r.h.s. of (|110|l liave a simple physical meaning: 
the first is the hard, medium-independent radiation H1U4|) reduced by the 
probability wi that one interaction of the projectile occurs in the medium. 
The second term describes the hard radiation component which rescatters 
once in the medium. The third term is the medium-induced Gunion-Bertsch 
contribution associated with the rescattering. For the case iV = 1, the 
above discussion that eq. (|106|l pattern interpolates between simple and 
physically intuitive limiting cases. In general, the A^-th order term in the 
opacity expansion is a convolution of the radiation associated to A^-fold 
scattering and a readjustment of t he p robabilities that rescattering occurs 
with less than N scattering centers'"^. 

4.2.2. Qualitative estimates of vs. quantitative calculations 

In general, the radiation spectrum lu^ displays an interference pattern 
between medium-induced and hard (vacuum) radiation amplitudes. For- 
mation time effects determine this interplay and allow for some qualitative 
estimates. In the simplest case, consider a hard partonic projectile which 
picks up a single transverse momentum fi by interacting with a single hard 
scatterer. An additional gluon of energy uj decoheres from the projectile 
wave function if its typical formation time fcoh = is smaller than the 
typical distance L between the production point of the parton and the 
position of the scatterer. The relevant phase is 




tcoh W 



(111) 



Thus one expects that the radiation of gluons is suppressed if their energy 
Lu is larger than the characteristic gluon energy 
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From this one can estimate the gluon energy spectrum per unit path length 
in terms of the coherence time fcoh, 

^ , , - T — ~as — . (113 

duj dz tcoh ^ 

This estimate can be compared to the fuh iV = 1 opacity term (|106|l inte- 
grated over transverse momentum up to |k| < x^- Introducing the kine- 
matic constraint in transverse momentum phase space, 

R = X^^cL, (114) 

the gluon radiation spectrum takes the forrr^l 

dl^=^ a^GR , , f°° r- sin{r) 

10— = 2 ("0^)7 / 5 

TT Jo ^ 

/ 1 1 A 



- = . (115) 
\r + l y/{{R/2-i) + r + 7)2 - ) 

The expression (jll5ll is plotted in Fig. El It shows qualitatively the same 
shape as the spectrum in the multiple soft scattering case. The position of 
the maximum of ^ '^-^(^-i) changes cx -i=. This is in accordance with the 
assumption that gluon radiation is suppressed for energies below which the 
characteristic angle of the gluon emission is of order one 



R Vr 

In the limit i? — s- oo in which the phase space co nstr aint is removed, 
one recovers some pocket formulas from the literature ESI. First, this limit 
satisfies the characteristic l/w-energy dependence of the estimate (I11H|I for 
sufficiently large gluon energies lo > Cuc, 

lim u! — 2 (no L) 7 / dr ■ 



^2^inoL){'^^P)^ ^^>"(117) 
I 4 77 for Wc < w 

Moreover, this expression allows to estimate the average parton energy loss 
for a single hard scattering. Re marka bly, this seems dominated by contri- 
butions from the region uj > (D^ p^ l ^^ l 

dujuj—, ~ ^^(noL)u^, log[i;M] , 

duj 2 

(118) 
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which appears to be logarithmically enhanced in comparison to the region 
LU < UJc for which 

lim / acuu! — ~ [noLjuJc- (119) 

R^oc Jq dcj TT 

However, the logarithmic enhancement of H118() does not persist upon closer 
inspection: All calculations so far worked for fixed coupling constant which 
may be justified if all momentum transfers are of the same order. In calcu- 
lating H118|l . however, the k-integration was extended up to the total energy 
E. If one takes the running of the coupling into account, this logarithmic 
enhancement factor is reduced to a log log factor. 
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4.2.3. Parameters in the opacity expansion 

In the multiple scattering approximation, the gluon energy distribution (|69|) 
depends on two quantities: the transport coefficient q and the in-medium 
path length L (we ignore the angular dependence in the following discus- 
sion) . In the opacity expansion, there are three parameters instead: the in- 
medium path length L, the transverse momentum squared /i^ transferred 
on average by a single scattering, and the inverse mean free path uq = j. 
The multiple scattering approximation describes the average momentum 
transfer per unit path length while the opacity expansion specifies in ad- 
dition by how many scattering centers this momentum is transferred on 
average. In the limit of multiple soft scattering, this additional information 
is redundant since 

fj.^ no L = q L , for Brownian motion. (120) 

However, deviations from Brownian motion arise due to the high trans- 
verse momentum tails of the elastic scattering cross sections |a(q)P = 
(/i'^(27r)^) / (7r(q^ -|- /i^)^). These tails lead to a logarithmically enhance- 
ment factor in the transport coefficient 

qL^noL ^ \a{q)\^ i cos^ ^ ~ {noL) fi^ \n M . (121) 



This allows to relate via ujc — (no-^)^cln y ^ the characteristic gluon en- 
ergies (|85|l and 1)112(1 which set the energy scale in the multiple soft and 
single hard scattering limits. For realistic values [/i > Aqcd and Q < E say] , 

Ini/^ <C 10, one concludes from Figs. ^ and Elthat the medium-induced 

V A* 

gluon energy distribution is significantly har der in the opacity approxima- 
tion than in the multiple soft scattering limiJ^. 



5. Applications 

Irrespective of the number of additionally radiated gluons, what matters for 
the medium modification of hadronic observables is how much additional 
energy Ai? is radiated off a hard parton. In this section, we first discuss 
the so called quenching weight which is the probability distribution P{AE) 
of the additional medium- induced energy loss. For independent gluon emis- 
sion, this probability is the normalized sum of the emission probabilities for 
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an arbitrary number of n gluons which carry away a total energy 



OO ^ 71 p 

n^^) - E n / 

n=0 ' U=l 



duj 



(122) 



Then we discuss how this probabihty can be used to calculate the medium 
modification of hadronic observables. 



5.1. Properties of Quenching Weights 

In gi 
partf 



In general, the quenching weight H122f) has a discrete and a continuous 

EI] 



P{AE) ^pa6{AE) +p{AE) . (123) 

The discrete weight po emerges as a consequence of a finite mean free path. 
It determines the probability that no additional gluon is emitted due to 
in-medium scattering and hence no medium-induced energy loss occurs. 

In order to determine the discrete and continuous part of H123f) . it is 
convenient to rewrite eq. 1)1221) as a Laplace transformation 121 



PiAE)^ j^±-V{v)e'^^''., 



V{v) = cxp 



oo 







, dl(u;) , 
duo—-^ (l 
duj 



(124) 



(125) 



Here, the contour C runs along the imaginary axis with Ret^ = 0. 

For the further discussion, it is useful to treat the medium-induced gluon 
energy distribution in eq. (|69|l explicitly as the medium modification 
of a "vacuum" distribution lS3 

rf/(tot) Jj(vac) 
OJ— ^LO ; +L0—. (126) 

duj aui dio 

From the Laplace transform (|124|l . one finds the total probability 

/•oo 

pi^°^){AE)^ dEP{AE- E)P^^'"'\E). (127) 
Jo 

This probability p(*°*) [AE) is normalized to unity and it is positive definite. 
In contrast, the medium-induced modification of this probability, P{AE), 
is a generalized probability. It can take negative values for some range in 
AE^ as long as its normalization is unity, 

/>00 /"OO 

/ dEP{E)=po+ dEp{E) = l. (128) 
Jo Jo 
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We now discuss separately the properties of the discrete contribution pq and 
the continuous one p{E). A CPU-inexpensive Fortran routinJ^is available 
for the calculation of these quenching weights. 
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Fig. 7. The discrete part po of the quenching weig ht c alculated in the multiple soft 
scattering limit as a function of R. Figure taken fromE^l 



5.1.1. Discrete part of the quenching weight 

The discrete part of the quenching weight is the n = term of eq. (|122ll . It 
can be expressed in terms of the total gluon multiplicity, 



Po = lim V{v) = exp [-N{uj = 0)] , 



(129) 



where the multiplicity N{ijj) of gluons with energy larger than uj emerges 
by partially integrating the exponent of H125|) . 



N{lo) 



duj' 



(130) 
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For the hmiting case of infinite in-medium path length, the total multi- 
plicity N{u;) diverges and the discrete part vanishes. In general, however, 
Po is finite. A typical dependence of po on model parameters is shown in 
Fig. for the radiation spectrum calculated in the multiple soft scattering 
limit. A qualitatively similar behavior is found in the opacity expansion. Re- 
markably, Po can exceed unity for some parameter range, since the medium 
modification uj-^ to the radiation spectrum H12t)|) can be negative. The 
value Po > 1 then compensates a predominantly negative continuous part 
p{AE) and satisfies the normalization (|128|l . This indicates a phase space 
region at very small transverse momentum, into which less gluons are emit- 
ted in the medium than in the vacuum. This effect is more pronounced for 
gluons than for quarks. 

5.1.2. Continuous part of the quenching weight 

The continuous part p{AE) of the probability distribution H123|l can be 
calculated numerically from the Mellin transform (|125|l . To facilitate the 
numerical calculation, one subtracts from P{v) the discrete contribution po 
which dominates the large-;/ behavior. 

Fig.lHlshows the continuous part of the quenching weight, calculated in 
the multiple soft scattering limit. In the opacity expansion, it looks qual- 
itatively similar. As expected from the normalization condition (|128|l . the 
continuous part p{AE) shows predominantly negative contributions for the 
parameter range for which the discrete weight po exceeds unity. With in- 
creasing density of scattering centers (i.e. increasing R — ^qL^) the proba- 
bility of loosing a significant energy fraction AE increases. The energy loss 
is larger for gluons which have a stronger coupling to the medium. This 
broadens the width of p{AE) for the gluonic case. 

In the multiple soft scattering approximation, an analytic estimate for 
the quenching weight can be obtained*' in the limit R ^ oo from the small- 
Lo approximation ujjr-ca -y= 



This reproduces roughly 1^ the shape of the probability distribution for 
large system size, but it has an unphysical large e-tail with infinite first 
moment J dee ^bdm™(£)- alternative analytic approach^ aims at fitting 
a two-parameter log-normal distribution to the numerical result for P{AE). 




where a = 




(131) 
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Fig. 8. The continuous part of the quenching weight 11231 . calculated in the multiple 
soft scatte ring limit for a hard quark (upper row) or hard gluon (lower row). Figure 
taken from 1221 



5.2. Quenching factors for hadronic spectra 

Assume that a hard parton looses an additional energy fraction IS.E while 
escaping the collision region. The medium-dependence of the corresponding 
inclusive transverse momentum spectra can be characterized in terms of the 
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quenching factor Q 



m 

^^^^^ da^^^{p^)ldp\ J V da^^^{p^)/dp\ J 

^JdAEPiAE){^-J^y . (132) 

Here, the last hne is obtained by assuming a power law fall-off of the p±- 
spectrum. The effective power n depends in general on p± . It is n ~ 7 for 
the pj^-range relevant for RHIC. Alternatively, instead of the quenching 
factor (|132|) . the medium modification of hadronic transverse momentum 
spectra is often characterized by a shift factor S{p±), 

da'^'^^jp^) da^-%pj_+S{p^)) 

dpi " d^i ' ^'''^ 

which is related to the shift S{p±) by 

0(p±)-exp|-^.5(pi)| . (134) 

Most importantly, since the hadronic spectrum shows a strong power law 
decrease, what matters for the suppression is not the average energy loss 
(A£') but the least energy loss with which a hard parton is likely to get 
away. One concludes that S{p±) < (AE) and depends on transverse mo- 
mentum ^. Fig. shows a calculation of the quenching factor p32|l in the 
multiple soft scattering limit. A qualitatively similar result is obtained in 
the opacity expansion. In general, quenching weights increase monotonically 
with p± since the medium-induced gluon radiation is independent of the 
total projectile energy for sufficiently high energies. At very low transverse 
momenta, the calculation based on (|69() is not reliable and the interpre- 
tation of the medium modification of hadronic spectra in nucleus-nucleus 
collisions will require additional input (e.g. modifications due to the Cronin 
effect). Fig. suggests, however, that hadronic spectra at transverse mo- 
menta p± > 10 GeV, can be suppressed significantly due to partonic final 
state rescattering. 

To quantify the sensitivity of the calculation to the low momentum 
region, Baier et al.l^ introduced a sharp cut-off on the R ^ oo gluon energy 
distribution which was varied between Wcut — and Wcut = 500 MeV. 
However, phase space constraints (i.e. finite R) deplete the gluon radiation 



46 



A. Kovner and U.A. Wiedemann 



CL 

^ 0.8 
0.6 

0.4 

0.2 

^0.8 
0.6 
0.4 
0.2 



q^=1 GeVVfm 
n=4, L=2fm n=4, L=5fm 



/ """ ^cut 


MS 

=500MeV 

= 300MeV ..:--^00^ 


Same cl 

R< c 

III 


jts with ^ — — ■ 



50 



100 

Pt(GeV) 



50 



100 

Pt(GeV) 



Fig. 9. The quenching factor il32i calculated in the multiple soft scattering limit. Upper 
row; calculation in the R — > oo-limit but with a variable sharp cut-off on the infrared 
part of the gluon energy distribution. Lower row: the same calculation is insensitive to 
infrared co ntrib utions if the finite kinematic constraint R = WcL < oo is included. Figure 
taken from EH 



spectrum in the soft region, see Fig. [3 As seen in Fig. IHl this decreases 
significantly the sensitivity of quenching factors to the uncontrolled infrared 
properties of the radiation spectrum. 
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5.3. Medium- modified fragmentation functions 

For an alternative calculation of the medium modification of hadronic spec- 
tra, one may determine the dependence of fragmentation functions on par- 
tonic energy loss. In general, hadronic cross sections are calculated by 
convoluting the parton distributions of the incoming projectiles with the 
product da^{z,Q'^) of a perturbatively calculable partonic cross section 
cr' and the fragmentation function Dj-^/ij{x,Q'^) of the produced parton, 
da^(z,Q2) ^ (^d^^ dyDh,q{x,Q^)dx. Here, x = Eh/ Eg, y = EJQ and 

z — Eh/Q denote fractions between the virtuality of the hard process Q, 
and the energies of the produced parton and resulting hadron. If the pro- 
duced parton loses with probability P{t) an additional fraction e = of 
its energy due to medium-induced radiation, then the hadronic cro ss sect ion 
is given in terms of the medium- modified fragmentation function I^EHI 



D^-;\x,Q^) = |'&P(,)^D,/,(^,Q2). (135) 



The hadronized remnants of the medium-induced soft radiation are ne- 
glected in the definition of (|135|) . However, these remnants are expected to 
be soft, and their inclusion would thus amount to an additional contribution 
to D^^l^\x, Q2) for X > 0.1 say. 

Fig. EH shows a calculation of the parton fragmentation functions 
_D^/^(x, Q^) from (|135(l using the quenching weights of Fig. IHland the LO 
KKP parametrization of Dfi/g{x, Q^). For this calculation, the virtual- 
ity Q of Dfi/q{x,Q^) is identified with the (transverse) initial energy Eg 
of the parton. This is justified since Eg and Q are of the same order, and 
Dh/q{x,Q^) has a weak logarithmic Q-dependence while medium-induced 
effects change as a function of e = « ^(^)- ^"-"^ ^ collision region 
expanding according to Bjorken scaling, t he tr ansport coefficient can be 
related to the initial gluon rapidity density 



1 dNS 
W 



That's what is done in Fig. ^] Interestingly, eq. (|136|l indicates how par- 
tonic energy loss changes with the particle multiplicity in nucleus-nucleus 
collisions. This allows to extrapolate parton energy loss effects from RHIC 
to LHC energies 
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Fig. 10. The LO KKP l^il fragmentation function u — > tt for no medium and the 
medium-modified fragmentation functions for different gluon rapidity densities (see eq. 
I13bi ') and L = 7 fm. Figure taken fromlSH 



In principle, the medium modified fragmentation function should be 
convoluted with the hard partonic cross section and parton distribution 
functions in order to determine the medium modified hadronic spectrum. 
For illustration, however, one may exploit that hadronic cross sections weigh 
D^f^^'^\x,Q'^) by the partonic cross section ci(T''/cip5_ ~ and thus 

effectively test x'^^^'P^^ d["!'"^\x, Q"^)^. The value n = 6 characterizes!^ 
the power law for typical values at RHIC {y/s — 200 GeV and p± ~ 10 
GeV). Thus, the position of the maximum Xmax of x^ d''j^°'^\x , Q^) corre- 
sponds to the most likely energy fraction Xma^Eq of the leading hadron. 
And the suppression around its maximum translates into a corresponding 
relative suppression of this contribution to the high-pj^ hadronic spectrum 
at p± ^ a^max-E'g- In general, the suppression of hadronic spectra extracted 
in this way is in rough agreements^ with calculations of the quenching 
factor (dSSl- 
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6. Appendix A: Eikonal calculations in the target light 
cone gauge. 

In section 121 of this review we have used the light cone gauge A~ = 0. In 
this gauge the gluon and quark distributions of the projectile wave function 
are simply expressible in terms of the gluon and quark number operators 
and we will therefore refer to it as the projectile light cone gauge (PLCG). 
The standard light cone gauge used in DIS calculations on the other hand 
is = 0, which facilitates simple expressions of the target distribution 
functions. In this appendix we show how the spectrum of emitted gluons 
is obtained in this standard light cone gauge, which we will refer to as the 
target light cone gauge (TLCG). 

Recall, that in PLCG the target is described as an ensemble of gluon 
fields with dominant component A^. The eikonal 5- matrix for the propa- 
gation of a charged parton through the target is given my the Wilson line 
eq.© 

VF(x,) = V exp{i J dz.T'^A+Cx,, z_)} (137) 

In the TLCG the yl+ component of the vector potential vanishes. Instead 
the chromoelectric fields in the target are given in terms of the transverse 
components Ai. The are obtained from by the gauge transformation 
from PLCG to TLGT^ZI 

A,(x„a:_) = iVHx„X-)d,V{x^,X-) (138) 

where 

F(x„x_) -Pexp{i / dz_T'^A+(x„z_)} (139) 

Note that the TLGT condition — does not fix the gauge unambigu- 
ously, but only up to residual gauge transformation which does not depend 
on x~ . The choice of the lower limit of the z_ integration in ea. (|139|l is 
equivalent to fixing this resi dual gauge freedom by imposing the condition 
djAj{'Ki,x- — oo) = qI^. With this choice we have 

F(x,,2;_ ^ -cx)) = 1, l/(x,,x_ -> oo) = H/(x,) . (140) 

Ea. l|138(l defines the vector potential Ai as two dimensional pure gauge, 
diA'^ - djAf - f^^A^A'^^ = 0. Moreover, at x_ ^ -l-oo the vector potential 
is genuinely (and not just two dimensionally) pure gauge. 
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Let us now consider scattering of a projectile with the wave function 
eq.(^ on the target described by ensemble of fields Ai of the form ea. (|138|l . 
Since the component of the vector potential vanishes, in eikonal ap- 
proximation the wave function of the projectile does not change while it 
propagates through the target. The outgoing wave function therefore is 
equal to the incoming one 

^'o«t = V'({ai,xJ) |{a„xj) . (141) 

{ai,Xi} 

This however does not mean that no scattering takes place. To calculate 
the scattering amplitude one has to project the outgoing wave function into 
the Hilbert space orthogonal to the wave function of the freely propagating 
system far "to the right" of the target, that is at x- — > +00. In the PLCG 
the target gauge field vanishes at both x_ — > — 00 and x- — > +00 and the 
freely propagating wave functions are identical "to the left" and "to the 
right" of the target. However in TLGT this is not the case. At x_ — > +00 
the target vector potential does not vanish, but is instead a pure gauge, eq. 
(|138|l . Therefore the freely propagating wave function at x^ — > +00 is not 
identical to that at a;_ — > —00, but is rather its gauge transform with the 
gauge transformation generated by the Wilson loop eq. H137|) . In particular 
the fields "to the right" of the target, A are related to the fields "to the 
left" of the target, a by 

A^W^aW + iW^dW . (142) 

Thus the free Fock basis "to the right" of the target is related to the free 
Fock basis "to the left" of the target by 

|{a,xJ)H = W^]^(x.)|{/3,x,})i . (143) 

The outgoing wave function ea. (|141|l is given in the basis |{Q!,Xi})L. On 
the other hand all observables at late time after scattering, including the 
number of emitted gluons must be calculated with respect to the basis 
\{a,Xi})B,- It is thus convenient to rewrite "if out using ea. (|143|l as 

*ont- VX{a^,xJ)nW^(xOc«,ft|{A,xJ)fl. (144) 

{ai,Xi} i 

However the same projectile freely propagating to the right of the target 
would have the wave function 

*/ree= ^ V^ll"* , X, }) | {a„ X, }) _R . (145) 

{Oi.Xi} 
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Thus we see that the outgoing wave function indeed differs from the 
freely propagating one, and thus the scattering is nontrivial and gluons are 
emitted in the final state. It is a somewhat curious feature of the TLCG 
that the nontrivial scattering amplitude appears entirely due to rotation 
of the free particle basis between early and late times Nevertheless it is 
obvious that the results of the calculation in this gauge are identical to 
those in PLCG as they should be. In fact from this point on all calculations 
are identical to those presented in section |21 as the interesting part of the 
outgoing wave function is given by eq.© with VE" free of ea. H145|l substituted 
for VPin and the "free" gluons are defined as states in the R Hilbert space, 
|{Q!,x})fl. 



7. Appendix B: Path integral formalism for the photon 
radiation spectrum 

In section rOl we discussed how to derive the non-abelian gluon radiation 
spectrum from the non-abelian Furry approximation H44(l . In this appendix, 
we present in more detail the derivation of the abelian analogue from the 
abelian Furry wave function H5()|l . The QE D photon radiation cross section 
in terms of Furry wave functions read*^ 

'^''^ "^"^ \MM\ (146) 



d{lnx) dp dk {2T:y 



ld^x.-\x,,.)...e-^^^,^-^.^ix,,.). (147) 



Here, e"'^'^' is the adiabatic switching off of the interaction term at large 
distances, discussed below eq. To simplify H147|) . we perform the fol- 

lowing steps: 

(1) Rotation of coordinate system Choose the momenta pi and p2 of the 
incoming and outgoing electron in the frame in which the longitudinal 
axis is taken along the photon: 



— 1 1 - X 
Pi = — k, P2=P k, (148) 

X X 

(2) z-dependent phase: The z-dependent phases of the Furry wave function 
l)5()|l combine in the radiation amplitude H147(l to an inverse photon 
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formation length 



q ^ pi - P2 - k = 



2{l-x)Ei' 

(3) Simplifying the spinor structure: The spinor structure 
Tr = ^/l — X u* (P2) D201 ■ ^Di w(pi) 



(149) 



(150) 



in the amplitude (|147|l can be simplified on the cross section level. The 
spin- and helicity-averaged combination F,*, take the simple form 



frV:,^[A-Ax + 2x^]^--^ +2mlx' 



(151) 



dv dr' 

(4) in-medium average: The cross section (|146|) contains products of the 
Green's functions H58() . These are averaged over the distribution of scat- 
tering centers in the medium, see (|52|) for the non-abelian case. These 
averages can be written in terms of the dipole cross section (|14|l : 



expjz J [C/(r(e),e)-C/(r'(0,e)] 



(152) 



For the products of Green's functions, this leads to 
{G{p,z';r,z\p2) G* {r' , z' ; p' , z\pi) ) = 



-P 1 
2'K i (z' — z) X 



X exp 



ip 



2{z' - z) 



X / Vrb exp <^ ^ 



where p = p{\ — x)x and 



{i-x){p-vf-{v'-p'y 



Pf[k) = P2- + Pi- 

z' — z z' — z 



(153) 



(154) 



After inserting the Furry wave function (|5U|I into the radiation amplitude 
H147|l . The radiation probability (|Myip) in terms of averages of Green's 



Gluon Radiation and Parton Energy Loss 



53 



functions reads 



oo 

(iMfip) = —^5- — Rc / dri dr dpdv2 dr[ dr' dp' dv'o I dz I dz' 

Ap{[l-xy J J J 

z 

X (G'(r2, 2;+; p, z'\p2) G* (r^, 2;+; r', z'\p2)) 
xf^r K'iGip, z'; r, z\p2) G*(r', z'; p\ z\pi)) 
x(G(r,z;ri,z_|pi)G*(p',z;r'i,z„|pi)) . (155) 

Equation H155(l is rearranged such that the photon emission in the amph- 
tude A^fi occurs prior to emission in the complex conjugate amplitude A^g. 
The opposite contribution is accounted for by taking twice the real part. Av- 
erages (...) of pairs of Green's functions effectively compare the paths r(^) 
and r'(^) of the electron in amplitude and complex conjugate am plitude , 
see eq. (|153ll . Using this average, one arrives at the final expressior 

d^a a„rn 2 



d(lna;) dk (27r)2 Ef (1 - x)^ 
X Re J dz J dz' exp 



X mz 



eiz + z' 



2(l-x) El 

Z- z 

J dri exp {— ik • ri} exp j' — y d^n{£,) a{xri) 

/C(z',r2 =0;z,ri|^), (156) 



4-4a; + 2a;2 d 



d 



? 2 
mi X 



Ax^ dri dv2 2 
where the abelian path-integral takes the form 

/C(z',rc(z');z,rc(2;)|M = £^i(l - x)x) 



Vvc exp ■ 



(157) 



The expression (|156|l provides a path-i ntegral formu lation of the Landau- 
Pomeranchuk-Migdal radiation spectrurrPSEnESEOI Graphically, as can be 
seen from eq. (I155|l . this radiation spectrum can be represented as 
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I 
1 
1 
I 




p ''2 


paf/i /n 


1 
I 
1 






[ll-— ' — ^ 

Pi=rrr; 


Pf(^)+ X r^l^) 


P2= V 

= p-r 


difference between 
the two paths 




P' 


r' 


path in 



Thus, for propagation from z_ = —00 up to the longitudmal photon emis- 
sion point z in Ma, the electron propagates with initial energy Ei in both 
amplitude and complex amplitude. Then it propagates with momentum 
E2 ~ Ei{l — x) in Ma but with Ei in TWg up to z', the photon emis- 
sion point in the complex amplitude, etc. The average H153() ensures that 
if propagation in Ma and A^g occur with the same energy, then the rela- 
tive distance between the paths r(^) and r'(^) does not change: physically, 
the radiation cross section counts only those changes in r{£^) — r'(^) which 
amount to phase shifts between amplitude and complex conjugate ampli- 
tude. Those are accumulated between the photon emission points z and z' 
in A^fi and M'^, respectively. 
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